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Cross-section measurements in the NOMAD experiment 
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The NOMAD experiment collected valuable neutrino data samples, matching both the large statistics of massive 
calorimeters and the reconstruction quality of bubble chambers. This paper describes the recent measurements 
of neutrino cross-sections on carbon target. The approach followed for cross-section modeling is also explained. 



1. Introduction 

The NOMAD experiment was designed to 
search for v T appearance from neutrino oscilla- 
tions in the CERN wide-band neutrino beam pro- 
duced by the 450 GeV proton synchrotron. The 
single-particle reconstruction and lepton identifi- 
cation capability of the NOMAD detector allowed 
the search for v T appearance in most of the lep- 
tonic and hadronic r decay channels !Tj and also 
to look for — > v e oscillations 2 . No evidence 
for oscillations was found. 

The high resolution and reconstruction quality 
makes the NOMAD data samples a valuable re- 
source for the neutrino physics community. The 
recent measurements of cross-sections and par- 
ticle production would help to clarify our un- 
derstanding of neutrino interactions with nuclei 
at intermediate energies. In addition, it must 
be noted the precise tracking in a light target 
provides a powerful tool to tune and validate 
Monte Carlo simulation programs for future ex- 
periments. This analysis activity can benefit from 
the beam and detector studies performed for the 
oscillation searches. 

2. Detector and data samples 

The NOMAD detector is described in detail 
in Ref. |3J- Inside a 0.4 T magnetic field there 
is an active target consisting of drift chambers 
(DC) |JJ with a fiducial mass of about 2.7 tons 



and a low average density (0.1 g/cm 3 ). The main 
target, 405 cm long and corresponding to about 
one radiation length, is followed by a transition 
radiation detector (TRD) for electron identi- 
fication, a preshower detector (PRS), and a high 
resolution lead-glass electromagnetic calorimeter 
(ECAL) [SJ. A hadron calorimeter (HCAL) and 
two stations of drift chambers for muon detec- 
tion are located just after the downstream part 
of the magnet coil. An iron-scintillator sam- 
pling calorimeter with a fiducial mass of about 
I7t (FCAL) is located upstream of the central 
part of the NOMAD target. The detector is de- 
signed to identify leptons and to measure muons, 
pions, electrons and photons with comparable res- 
olutions. Momenta are measured in the DC with 
a resolution: 



0.05 0.008 x p[GeV/c 



P 



where L is the track length and p is the momen- 
tum. The energy of electromagnetic showers, E, 
is measured in the ECAL with a resolution: 



-f = 0.0! 



0.032 
y/E[GeV] ' 
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The relative composition of Charged Current 
(CC) events in NOMAD is estimated to be 
v ll ,:v^'.v s ,:v e = 1.00: 0.0227: 0.0154: 0.0016, 
with average neutrino energies of 45.4, 40.8, 57.5, 
and 51.5 GeV, respectively. Neutrinos are pro- 
duced at an average distance of 625 m from the 
detector. 
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The NOMAD experiment collected data from 
1995 to 1998. Most of the running, for a total 
exposure of 5.1 x 10 19 protons on target (pot), 
was in neutrino mode. This resulted in three 
distinct data samples, according to the differ- 
ent targets: 1.3 x 10 6 z^ CC interactions from 
the drift chambers (mainly carbon), 1.5 x 10 6 ^ 
CC interactions from the region of the magnet 
coil (mainly aluminium) located in front of the 
DC and 1.2 x 10 7 ^ CC interactions from FCAL 
(iron). 

3. Modeling of inelastic cross-sections 

3.1. Structure functions at high Q 2 

At large momentum transfer structure func- 
tions are described as series in Q~ 2 on the basis of 
the operator product expansion of the correlator 
of the weak current (twist expansion) : 



F t (x,Q 2 



TMC 



F, 



Q 2 
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(i) 



where i = T, 2, 3 refers to the type of the struc- 
ture function, F i TMC are the leading twist (LT) 
terms corrected for the target mass effects and 
h[ are the higher twist (HT) terms of twist t. 
The target mass corrections are computed using 
the approach of Ref.[Hl (see also f° r the treat- 
ment of the thereshold problem at x — > 1). The 
calculation of the leading twist is performed in the 
NNLO approximation and we include additional 
phcnomenological terms up to twist-6. We use 
PDFs and HT based on Ref.|§], obtained from 
dedicated fits optimized at low Q 2 ^D] and in- 
cluding additional data from (anti)neutrinos CC 
(NOMAD, CHORUS 25, and NuTeV g2]) and 
charged lepton (JLab) Deep Inelastic Scattering 
(DIS), as well as from Drell-Yan production (E605 
and E886). 

3.2. Low Q 2 structure functions 

In the low-Q region (anti)neutrino cross sec- 
tions are dominated by the longitudinal struc- 
ture function Fl and the latter is driven by the 
axial-current interactions. The structure function 
Ft vanishes as Q 2 at low Q 2 . This behaviour 
is similar to the charged lepton case and holds 



for both the vector and the axial-vector contri- 
butions. However, in the longitudinal channel 
the low-Q behavior of the vector and axial- vector 
parts are different. 

The conservation of the vector current (CVC) 
suggests q^W^v — for the vector current part 
of hadronic tensor. From this condition we con- 
clude F^ vanishes faster than Fj, at low Q 2 and 
FlI^t ~ Q 2 - This behavior is similar to the 
charged-lepton case. 

In contrast to the vector current, the axial- 
vector current is not conserved. For low momen- 
tum transfer the divergence of the axial-vector 
current is proportional to the pion field (Partially 
Conserved Axial Current or PCAC) 



3A ± = um 2 y 



(2) 



where m„ is the pion mass and f n = 0.93™^ is 
the pion decay constant and ip^ is the pion field 
in the corresponding charge state. We introduce 
explicitely a PCAC contribution to F^\ 



F A =7 3 F PCAC /pcAc(g2) 



Ft 



(3) 



where 7 = (1 + 4x 2 M 2 /Q 2 ) 1 / 2 , Fl° AC = f 2 a^/ir 
and ov = (t^{s,Q 2 ) is the total cross section 
for the scattering of a virtual pion with mo- 
mentum q and the center-of-mass energy squared 
s = (p + q) 2 . The last term F A is similar to F^ 
and vanishes as Q 4 . Since the PCAC contribu- 
tion is expected to vanish at high Q 2 we introduce 
a form factor / PC ac(Q 2 ) = (1 + Q 2 /M 2 CAC )- 2 7 
where the dipole form is motivated by meson 
dominance arguments. It is important to note the 
pion pole does not directly contribute to structure 
functions and hence the mass scale controlling the 
PCAC mechanism, Mpcac> cannot be the pion 
mass itself, but is rather related to higher mass 
states like a±, pit etc. We use the simple assump- 
tion Mpcac = rn ai |T5|. 

The structure functions F T and F L = F^ + F^, 
which are vanishing for Q 2 — > like in the 
charged lepton case, are parameterized as smooth 
interpolations between the high Q 2 regime cal- 
culated from Equation ^ and the Q 2 — > pre- 
dictions derived from current conservation argu- 
ments |TD]|r3l. We choose the value Ql = 1 GeV 2 
as matching point for the twist expansion. In the 
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region < Q 2 < 1 GeV 2 we use cubic splines 
calculated for fixed x values. The coefficients of 
such functions are fully determined by the condi- 
tion both functions and derivatives should match 
with the twist expansion at Q\. Figure illus- 
trates the interpolation procedure for F% on pro- 
tons in charged lepton scattering. 
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Figure 1. Interpolation of structure functions in 
the region < Q 2 < 1 GeV 2 . The example given 
in the plot refers to F2 for charged lepton scat- 
tering on protons at x=0.07 (see text for details). 



The different Q 2 dependence of various terms 
in Equation ^ allows to disentangle higher twist 
contributions, which are parameterized as smooth 
functions (splines) of x. Our results indicate the 
twist-6 term oc 1/ Q 4 is important in order to de- 
scribe the ratio R of longitudinal to transverse 
cross-sections at low Q 2 . This is shown in Fig- 
ure |2] together with SLAC and JLab data. We 
derive the ratio R directly from the parameteri- 
zations of Fl and Ft ■ 

From the relation F 2 = (F L +F T )/j 2 and © it 
follows that the structure function F2 at low Q 2 
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Figure 2. Impact of higher twist terms on R = 
pl/pt for charged lepton scattering. Data from 
SLAC and JLab are shown for comparison. 



is dominated by the nonvanishing F^ CAC term 
(Figure |2J. It is important to note that since 
F T — > and F L — > F^ CAC in the limit of van- 
ishing Q 2 the ratio R — Fl/Ft is divergent for 
neutrino interactions. This is substantially dif- 
ferent from the scattering of charged leptons for 
which R is vanishing as Q 2 . 

The determination of LT and HT terms is 



performed from all available data with Q 2 > 
0.5 GeV 2 and W > 1.9 GeV. It is interesting 
to check the extrapolation of DIS structure func- 
tions into the resonance region. The results are 
consistent with the duality principle, as can be 
seen from Figure^ where the integral of the dif- 
ference between the recent JLab resonance data 
and the average DIS predictions is consistent with 
zero. 

3.3. Nuclear and electroweak corrections 

A detailed calculation of nuclear corrections to 
structure functions is performed JIJ^J. The 
model takes into account a number of different 
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Figure 3. The PC AC contribution to the neutrino 
structure function F2 at x = 10~ 4 for different 
targets. 



Figure 4. Comparison between the extrapolation 
of D1S F2 structure function on proton and the 
recent JLab data in the resonance region. 



effects including nuclear shadowing, Fermi mo- 
tion and binding, nuclear pion excess and off-shell 
correction to bound nucleon structure functions. 
The off-shell effect and the effective scattering 
amplitude describing nuclear shadowing are ex- 
pressed in terms of few parameters which are 
common to all nuclei and have a clear physical in- 
terpretation. The parameters are then extracted 
from statistical analysis of data from charged lep- 
ton scattering on a wide range of nuclear targets. 

The treatment of Fermi motion, binding, 
off-shell effect and nuclear pion correction in 
(anti)ncutrino interactions is similar to the one 
in charged lepton scattering. The main differ- 
ences are related to the impact of the axial- vector 
current on coherent nuclear effects and are more 
evident at low Q 2 [Tl). 

One-loop electroweak effects are taken into ac- 
count 53 as corrections to the parton distri- 
butions used in the structure function calcula- 
tion. The initial quark mass singularities from 
QED corrections are subtracted within the MS 



scheme since they are effectively incorporated 
in the quark density functions. Results are 
cross-checked with a second independent calcu- 
lation 

3.4. Neutrino fluxes 

Two independent approaches are used to cal- 
culate neutrino fluxes. The first method is based 
on a simulation of the West Area Neutrino Fa- 
cility 7 . The calculation of particle production 
rates from the interaction of primary protons on 
Be target is performed with a recent version of 
FLUKA fI2| , further modified to take into ac- 
count the cross-sections measured by the SPY 
and NA20 experiments. These particles are then 
propagated through the beam line taking into ac- 
count the material and magnetic fields they tra- 
verse. Predictions are then validated by compar- 
isons with NOMAD data. 

A second method relies directly upon the CC 
events reconstructed in the NOMAD detector. In 
particular, events with low energy associated to 
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the hadronic system (y < 3 GeV) are analyzed to 
this purpose since the corresponding differential 
distribution dN / dv is proportional to the flux up 
to a small correction factor (~ 10%). 

3.5. Monte Carlo simulations 

The Monte Carlo (MC) simulation of neutrino 
DIS interactions is based on LEPTO 6.1 ^H] an d 
JETSET 7.4 Q1] packages, followed by a full 
GEANT3 [23] propagation to model the detector 
response. The neutrino cross-sections are param- 
eterized according to the model described in the 
previous Sections. The simulation of resonance 
production is performed according to Ref. j21| . 

We do not include the parton shower treatment 
from JETSET. The reinteractions of hadrons 
with surrounding nucleons in target nuclei are 
described within the DPMJET |22 package. A 
detailed tuning of fragmentation parameters has 
been performed with hadrons reconstructed in 
CC interactions |23| . 

4. Measurement of inelastic cross-section 

The present knowledge of neutrino cross- 
sections is rather nonuniform. In the region 
E v > 30 GeV, where data from the large massive 
calorimeters (CCFR, NuTeV) are available, the 
uncertainty is about 2%. This increases to about 
20% at lower energies, due to the limited statis- 
tics of bubble chamber experiments. Measure- 
ments of both the total (Jqq and the differential 
da^ c / 'dxdy cross-sections for i>a on carbon are 
performed in NOMAD. One of the primary goals 
is to constrain the (anti)neutrino cross-section 
model and to reduce the corresponding system- 
atic uncertainties. This also provides the first 
measurement of inclusive cross-section on carbon 
target, with (Q 2 ) ~ 13 GeV 2 . 

The absolute normalization is obtained from 
the world average cross-section on isoscalar tar- 
get in the energy range 40 < E v < 300 GeV. The 
measurement is performed in the kinematic re- 
gion v > 3 GeV and Q 2 > 1 GeV 2 . The analysis 
of the lower Q 2 region down to about 0.3 GeV 2 
is currently being finalized. 

A comparison of the measured differential 
cross-section on carbon with the model predic- 



tions (Sec. El is illustrated in Figure0for E = 85 
GeV, indicating a good agreement in the entire 
kinematic region. 

The NOMAD results are also consistent with 
the measurements performed by CHORUS and 
NuTeV on different targets [T3|. Only in the re- 
cent NuTeV data a slight excess is observed for 
x > 0.5 H. The low Q 2 data from CHORUS 
support our treatment of the PCAC contribution 
to neutrino structure functions. 
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Figure 5. Differential neutrino CC cross-section 
on carbon at E — 85 GeV. The full circles show 
the NOMAD measurement while the curve with 
open circles is the prediction from our cross- 
section model. 
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5. Measurement of quasi-elastic cross- 
section 

The reconstruction and identification of the re- 
coiling proton track allows a measurement of the 
quasi-elastic cross-section v^n — > pTp on carbon 
in NOMAD. A kinematic analysis based upon a 
three-dimensional likelihood function is used to 
reject backgrounds from DIS and resonance pro- 
duction. Overall, about 8000 two track events are 
selected for the cross-section determination in the 
energy range 3 < E v < 100 GeV. A complemen- 
tary sample of about 16000 single track events is 
used as a cross-check of the reconstruction effi- 
ciency. 

In order to reduce systematic uncertainties, 
NOMAD is measuring the ratio of quasi-elastic 
cross-section with respect to DIS (W 2 > 4 GeV 2 ) 
processes. The absolute normalization is pro- 
vided also in this case by the world average cross- 
section on isoscalar target. 

The simulation of the quasi-elastic events 
is based upon the formalism by Llewellyn- 
Smith |26| and the axial form factor is param- 
eterized in the conventional dipole form. Nuclear 
effects and the Pauli suppression factor are imple- 
mented following a simple Fermi gas model. Final 
state interactions are described by the DPMJET 
package The main systematic uncertainties 
are related to the understanding of nuclear cor- 
rections and are currently being finalized. 
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